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We investigate the magnetic response of normal-metal-superconductor proximity systems for 
arbitrary concentrations of impurities and at arbitrary temperatures. Using the quasiclassical theory 
of superconductivity a general linear-response formula is derived which yields a non-local current- 
field relation in terms of the zero-field Green's functions. Various regimes between clean-limit 
and dirty-limit response are investigated by analytical methods and by solving the general formula 
numerically. In the ballistic regime, a finite mean free path reduces the non-locality and leads to a 
stronger screening than in the clean limit even for a mean free path much larger than the system size. 
Additionally, the range of the kernel describing the non-locality is strongly temperature-dependent 
in this case. In the diffusive limit we find a cross-over from local to non-local screening, which 
restricts the applicability of the dirty-limit theory. 



I. INTRODUCTION 

A normal metal in good metallic contact to a super- 
conductor acquires induced superconducting properties. 
The basic features of this proximity effect were already 
well understood in the sixties. One of these properties is 
the diamagnetic screening of an applied magnetic Jield, 
which has been studied in a series of experimentaluTu and 
theoretical worksH~@. 

Already early in this development it was recognized 
that the relevant length scale governing the supercon- 
ducting correlations is given by the thermal and impu- 
rity dependent coherence lengths in the normal metal. 
The thermal length is given by £t = /2ttT in the clean 
limit (mean free path I — > oo) and £d = m the 

dirty limit (I — ► 0). The finite thickness d of the normal- 
metal layer is an intrinsic geometric length scale of the 
proximity effect. The interplay between the three length 
scales Z, and d is relevant for the behavior of the mi- 
croscopic quantities such as the spatial decay of the pair 
amplitude or the spectral density of states. In this pa- 
per we investigate the range between the clean and the 
dirty limiting cases and find several intermediate regimes 
of interest, differing by the relative magnitudes of £r, I, 
and d. 

The theory of linear diamagnetic response of clean and 
dirty N-S graximity systems has already been studied 
extensivelyEMtj. While the dirty-limit theory was-ipund 
to be in agreement with early experimentaLworkli3, the 
samples studied in more recent experimentsErEj fail to be 
described satisfactorily by either the clean or the dirty 
limit. They exhibit the qualitatively different|-|behavior 
of an intermediate regime. In a previous papem we were 
able to fit an experiment in the low-temperature regime 
with the dirty-limit theory. However, the data for tem- 
peratures T > v-p/d could not be reproduced by the 



dirty limit theory. These experiments also show clear 
deviations from the clean-limit theory, even ifrfhe finite 
transparency of the interface is accounted forEm. On the 
other hand, the breakdown field seen in the nonlinear 
response to a magnetic field of some of the same sam- 
ples was found to agree fairly well with the clean-limit 
theoryul. In this paper, guided by a numerical study, we 
classify the intermediate regimes and show how the qual- 
itative discrepancies between theory and experiment can 
be resolved. In recent experiments on relatively rClffl-P 
samples a low-temperature anomaly was reportecOtm. 
The nearly perfect screening at T w vp/d was found to 
be reduced as the temperature was lowered further. This 
re-entrance effect has not been explained up to now. We 
will not address this problem directly here, but rather 
provide an understanding of those facets of the proximity 
effect which are the necessary basis for further investiga- 
tions. 

From a theoretical point of view there is a major qual- 
itative difference between the magnetic response in the 
dirty versus the clean limit. In the dirty limit the current- 
field relation is local and the screening can be almost 
complete. This is in strong contrast to the clean limit, 
where the current-field relation is completely nonlocal 
and the current depends on the vector .potential inte- 
grated over the whole normal-metal layem. As a conse- 
quence, there is an overscreening effect, i. e., the mag- 
netic field reverses its sign inside the normal metal, and 
the magnetic susceptibility is limited to 3/4 of that of 
a perfect diamagnet. In the intermediate regime, the 
superfluid density and the range of the current-field rela- 
tions, both diminishing with decreasing mean free path, 
are shown to affect the screening ability in a contrary 
way and thus compete in the magnetic response. 

We investigate the magnetic response at arbitrary im- 
purity concentrations starting from the quasi-classical 
Green's functions in absence of the fields, which are dis- 
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cussed in Section ||. On this basis, we develop a theory 
of linear current response in Section [II. We produce a 
general result with the well-known structure: the current 
density j(x) is given by a convolution of a kernel K(x, x') 
with the vector potential A{x). The kernel is given ex- 
plicitly in terms of the Green's functions in the absence 
of the fields. Our formula easily yields the basic consti- 
tutive relations of the LondonEJ- and PipparcO-type for 
a superconductor. With the help of this kernel we calcu- 
late the magnetic susceptibility of a proximity system at 
arbitrary impurity concentrations in Section [V. We find 



# 2 0f, x) + f(v F , x)f\v Fl x) = 1 



(2) 



that the impurities have non-trivial consequences on the 
magnetic response. The range of the integral kernel can 
be strongly temperature dependent, and can be given by 
either £t or I. In particular, we show that even for I 
considerably larger that d the spatial dependence of the 
integral kernel strongly enhances the magnetic response, 
as compared to the clean limit . 



II. QUASICLASSICAL EQUATIONS AND 
PROXIMITY EFFECT 

The basic set of equations appropriate for describing 
spatially inhoraogeneous superconductors was d&weloped 
by EilenbergertJ and by Larkin and Ovchinnikovt^l (for a 
recent collection of papers on the quasiclassical method, 
seellll). They are transport-like equations for the qua- 
siclassical Green's functions, i.e., the energy- integrated 
Gorkov Green's functions, that are derived from the 
Gorkov equations under the assumption that the length 
scales relevant for superconductivity arc much larger 
than atomic length scales. We treat the presence of elas- 
tic impurities within the Born approximation (the full 
T-matrix-fopasialism has been shown to lead to quantita- 
tive changeal3) . The Eilenberger equations take the form 
(e=N) 



-v F (V + 2ieA(x))f(v F ,x) 



(1) 



(2u + -(g(x)))f(v F , x) - (2A(x) + -(f(x)))g(v F ,x) 

T T 

v F (V - 2ieA)p(v F ,x) = 

{2u + -(g(x)))p(v F ,x) - (2A(x) + -(p (x)))g(v F , x) 

T T 

-v F Vg{v F ,x) = 

(A(x) + ±-{f(x)))(f(v F ,x)-p(v F ,x)) 

These are three coupled differential equations for the nor- 
mal (diagonal) Green's function g and the anomalous 
(off-diagonal) Green's functions / and /t. They depend 
on Matsubara frequency ui — TrT(2n+l), the elastic scat- 
tering time r = l/v F , and the Fermi velocity «p, (...) 
denoting the average over the Fermi surface (h = c = 1 
throughout). The superconducting order parameter A is 
taken to be real. We note that the w-dependence of the 
Green's functions has been omitted in our notation. The 
Green's functions obey the normalization condition 



and fulfill the symmetry relations g*(—v F , x) = g(v F , x) 
and f*(~v F ,x) — f^(v F ,x) . The current is given by 



i2ep F m sr^ i 
j{x) = T > (v F g(v F ,x)} 

7T 1 — ' 

ul>0 



(3) 



and depends only on the imaginary part of g, due to the 
above symmetry relations. 

In this paper, we consider a system shown in Fig. |l| 
consisting of a normal-metal layer of thickness d, which 
is in ideal contact with a semi-infinite superconductor. 
A magnetic field (0, 0, H) is applied parallel to the metal 
surface, producing screening currents (0,j(a;),0) along 
the surface, which depend on the coordinate x. The pair 
potential is taken to be a step function A(x) = A6(-i). 
Assuming a thickness d 3> £o = v f /2itT c , we can ne- 
glect the self-consistency of the pair potential. Further- 
more, we assume specular reflection at the normal-metal- 
vacuum boundary. 
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FIG. 1. Geometry of the proximity model system. The 
thickness of the superconductor is assumed to be much greater 
than £o, the pair potential is taken real and assumed to fol- 
low a step function: A(x) = A9(-i). In our gauge, the 
screening current and the vector potential are parallel to the 
NS interface. The interface is assumed to be perfect and the 
normal-metal-vacuum boundary to be specularly reflecting. 



In the absence of external fields (we denote the cor- 
responding Green's functions by go, fo and Jo) Eqs. (|l|) 
reduce to 



- v x — f (vx,x) = 2uj(x)f Q (v x ,x) -2A(x)g (v x ,x) (4) 
dx 



v x — fo(v*,x) = 2Q(x)fl(v^x) -2A(x)g (v x ,x) . 



We have introduced the effective frequency uj(x) = lu + 
(go(x))/2r and pair potential A (x) = A(x) + (fo(x)) /2t. 
Equations (||) imply that fo(v^,x) = f^(—v x ,x) and, 
since (fo(x)) = (fo(x))*, for real A we obtain a real 
(fo(x)), too. 
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Depending on the relative size of the thermal length 
£ti the mean free path I, and the thickness d we distin- 
guish the ballistic, the dirty and the intermediate diffu- 
sive regime that are discussed in the following subsec- 
tions. These regimes are also shown in Fig. 0. 




FIG. 2. Dependence of the magnetic response on thermal 
length £t = vf/2-kT and mean free path I. In the ballis- 
tic regime I S> min{^T,d} we distinguish three regions: (a) 
the clean limit with infinite range of the kernel exhibiting a 
reduced diamagnetism (overscreening), (b) the quasi-ballistic 
limit with finite range £t increasing screening at large tem- 
peratures, (c) the ballistic limit where the finite range I en- 
hances the screening although / 3> d. In the diffusive regime 
I <C ^T,d the range of the kernel is given by I. The dirty 
limit with nearly isotropic Green's functions is restricted to 
I <S £o,£t,g(. Note that the current-field relations can still be 
local or non-local depending on the relative size of penetra- 
tion depth and mean free path. For comparison, the conven- 
tionally assumed border line between clean and dirty limits 
{I = £t) is indicated by a dotted line. 



A. Ballistic regime 



At temperatures above the Andreev temperature, T ^ 
T4, only the first Matsubara frequency u> — nT is rel- 
evant and the decay of the /-function is governed by 
£t = vp/2irT. Ta determines the temperature at which 
the /-function acquires a finite value at the outer bound- 
ary. 

An estimation using the Eilenberger equation Eq. (^) 
easily shows that the clean-limit solution is valid for 



/ > dexp(2d/£ T ) if £r < d 
I > d if £ T > d 



(6) 
(7) 



We note that this includes the region d <C I <C £t, the fi- 
nite thickness preventing the small mean free path I <C £t 
of becoming effective. In the remaining part of the ballis- 
tic regime, see Fig. ||, the full solution is not known, but 
we may produce an approximate solution, which charac- 
terizes well the numerical results found below. Limiting 
ourselves to £t <C d allows us to consider the Green's 
function for the first Matsubara frequency u> — nT only. 
We restrict ourselves to the forward direction v x = +vp. 
From Eq. (Q) we find that /o ~ 2exp(— x/£t) remains 
unchanged as in (||), and / '<1 and 1 — g <C 1 obey 
the approximate equations, 



d 

dx 
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/o f W = -y(/oW) 



(8) 



I* 



9o(x)) 



•j(fo(x))f(x) 



with the approximate solutions 



1 - 9o{x) 
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(9) 
(10) 



Here we have used (/o) ~ f/2, which is valid since 
/o(— I f x |) -C fo{\v*\)- Interestingly, while the induced 
superconducting correlations as described by (fo(x)) re- 
main unchanged as compared to the clean limit (|5|), the 
values of /q and 1 — go are of order £t/1 rather than ex- 
ponentially suppressed as in (^). This is of importance 
for the current response as we show below. 



The ballistic regime is limited by I min{^T,d}, 
which ensures a ballistic propagation of the electrons 
over the thickness or the thermal length of the nor- 
mal layer, respectively. As a limiting case, for I — ► 00 
(clean limit) Eq. (ji|) may be solved analytically. For 
T c 3> Ta = Vfi/2iTd, the solution in the normal metal 
takes the formQ 



Md-x) 



fo{v*,x) = e 

fo( v x,x) = f(-v X) x), Xd = 



I cosh Xd 

2ud 



(5) 



B. Dirty limit 



If impurity scattering dominates, as described by 

(90} / T ^ w aim (/oVfn^ ^ < 1- @ can ^ e re duced 
to the Usadel equationtj for the isotropic part (f Q (x)). 
Assuming uj A the solution in the normal metal takes 
the form 



(/o0z)> =co8h(J^(d- a ;))/cosh(J^d) , (11) 
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where D = uf.r/3 is the diffusion constant. Equation 
( |TT| ) shows that the important energy scale is the Thou- 
less energy — D /2-Kd 2 . The coherence length in this 
case is £d(T) = (L>/27rT) 1 / 2 , which reflects the diffusive 
nature of the electron motion. 

In the normal metal I -C £t, d are necessary conditions 
for the Usadel theory to be valid. However, as the numer- 
ical results will confirm below, the Usadel theory in the 
normal metal may not be applied without considering 
the superconductor inducing the proximity effect. The 
application of the Usadel equations requires the Green's 
functions to be nearly isotropic, which in the supercon- 
ductor is only fulfilled for (/o)/ T 3> A. The validity of 
the Usadel theory (in the absence of fields) is thus re- 
stricted to I <C d, £r, and I <C £o = vf/2ttT c , the dirty 
limit, see Fig. ^. 



C. Intermediate diffusive regime 



Now we relax all restrictions on the mean free path 
and investigate the regime between the ballistic regime 
and the dirty limit. Eq. (||) can be formally decoupled 
using the Schopohl-Maki transformationE] 
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FIG. 3. Spatial dependence of the anomalous Green's 
function fo(x,v x ) in a proximity normal-metal layer. The 
thickness is d = 10vf/T c and the frequency is to — Vp/d. The 
mean free path is I = d in (a) and I = O.ld in (b). The 
insets show the angular dependence /(vf x ) at the positions 
indicated by the arrows. We note that in the Usadel theory 
the angular dependence would be given by a linear function. 



fo(Vx,x) 



ao(vK,x) = — j- -,al{v x ,x) = — f -, (12) 



leading to the Riccati differential equations 



v x — a (v x ,x) = 2u}(x)a (vx, x) + A(x)(ag(w x , x) - 1) 



J 2 i 



Eqs. ( |13t ) provide the basis for a (stable) numerical so- 
lution. We have determined the impurity self-energies 
self-consistently by an iteration procedure starting from 
the dirty-limit expression. Representative results of the 
numerical calculation are shown in Fig. ||. 



We have chosen the frequency u> = VF/d and a mean 
free path oil = din Fig. |(a), / = O.ld in Fig. |(b). Note 
the distinction of the /-function for forward (i>f x = Vp, 
solid line) and backward propagation (wp x = — ^f, dashed 
line). As we cross over from the ballistic to the diffu- 
sive regime, the backward propagating branch changes 
from a monotonically increasing /-function of x for I = d 
to a decaying /-function for I = O.ld. In the ballistic 
case the backward moving electrons carry superconduct- 
ing correlations only after reflection from the normal- 
metal-vacuum boundary x = d. In the diffusive case, 
(13l)he backward propagating /-function is generated by the 
impurity scattering from the forward branch, thus tak- 
ing the same functional dependence of x, see Fig. § This 
behavior is illustrated in the insets, where / is plotted as 
a function of v x for fixed positions. The sharp features 
present in Fig. ||(a) are washed out by impurity scatter- 
ing in^(b). Remarkably, however, they are still far from 
the dirty- limit behavior, for which /o(t'Fx) in the insets is 
expected to be a straight line. According to conventional 
wisdom I = O.lfF /lu <C £t would indicate the dirty limit. 
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Considering that we have chosen d = 1Qvf/T c , implying 
I ss v-p/A, we notice that the dirty-limit condition is not 
fulfilled in the superconductor. The anisotropy of the 
/-function present in the superconductor by proximity 
induces an anisotropy inside the normal metal, which is 
not accurately described by the Usadel theory. For the 
dirty-limit theory to be valid in the normal metal, the 
superconductor has to be dirty as well. 

We note that in our present calculation we have as- 
sumed the same mean free path in the superconductor 
and the normal metal. Allowing for different mean free 
paths would affect the definition of the dirty and the in- 
termediate regime, and we do not further investigate this 
question here. 



III. LINEAR-RESPONSE KERNEL 

In this section we derive the general linear-response 
kernel (^|]) of a normal-metal-superconductor sandwich 
in terms of the Green's functions in absence of the fields. 
We consider the quasi-one-dimensional system shown in 
Fig. [l], assuming a superconductor of thickness d s and a 
normal metal of thickness d. The magnetic field is applied 
in z-direction as described by the gauge A — A(x)e y . To 
calculate the linear diamagnetic response, we separate 
the Green's functions into its real and imaginary parts, 
where the imaginary part is of first order in A and the 
real part of zeroth order: 

f(v F ,x) = f (v x ,x) + if 1 (v x ,v y ,x) 

P(v F ,x) = f^(v x ,x)+ift(v x ,v y ,x) (14) 

g(v F ,x) = ga{v x ,x) + igi(v x> v y , x) . 

The zeroth-order parts obey Eq. (Q) discussed in the pre- 
vious section. The first-order parts of Eq. (|l|) read 

-v x -^-fi(v x ,v y ,x) = 2w(x)f 1 (v x ,v Y ,x) (15) 
-2A(x)gi(v x ,v y ,x) + 2ev y A(x)f (v x ,x) 

Vx-jj^fl(Vx,Vy,X) = 2ul(x)fl(Vy i ,Vy,x) 

-2A(x) gi (v x , v y , x) + 2ev y A(x)fl(v x , x) . 
where u>(x), A(x) were given after Eq. (@) and 

gi(v x ,v y ,x) = (16) 
fojvx, x)ft{v x , v y ,x) + h(v x , v y , x)fl(v x , x) 
2^0 («x>£) 

follows from the normalization (^). We now apply the 
Maki-Schopohl transformation defined in ( |l^ ) to the full 
equations of motions ([!]) . After linearization we obtain 



fl(V X ,Vy,X) = 2 

fl(v x ,v y ,x) = 2 



ai(v x ,v y ,x) - a 2 ) (v x ,x)a\(v x ,v y ,x) 

(1 + a (v x ,x)al(v x ,x)) 2 
a\(v x ,Vy,x) - al 2 (v x ,x)a 1 (v x ,v y ,x) 



(1 + a (v x ,x)al(v x ,x)) 2 
and Eqs. (fla) are decoupled into 



v x d 

--- ai (v x ,v y ,x) 



(17) 



Q(x) + A(x)a (v x , x) ai(v x , v y , x) + ev y A(x)a (v x , x) 

Vx_d_ 
2 dx 

Q(x) + A(x)a) (v x , x) a\(v x , v y , x) + ev y A(x)a\ ) {v x , x). 



a\{v x ,Vy,x) = 



For a more general form of these equations that has 
been used to treat the linear electromagnetic response 
of vortices numerically, see ||l]. We will now proceed 
analytically. As a consequence of Eq. (]l7]) we find 
Oi(— Vy) = —ax(vy) and the same for a\, which leads 
to (fi) = (ft) = 0, as was already noted above. Further- 
more, since a\(v x ) — cii(—v x ), we only have to consider 
one of the two equations (e.g., the first one). Equation 
( |l7j ) is an inhomogeneous first-order differential equation, 
which can be integrated analytically. Assuming that / 
and p do not change sign with the help of (0) the solu- 
tion can be written as 



ai(v x ,v y ,x) = c(v x ,v Y ) 
2ev v 



m(v x ,x,x ) 



v x fl(v x ,x) 
where 



1 



fl(v x ,x) 
g (v x , x')]m(v x , x, x')A(x')dx' 



(18) 



(2 f A(x") 

i(v x , x, x ) = exp — / —7 

y \vJ x fl{v x ,x" 



') 



■dx" 



(19) 



In this equation xq is an arbitrary reference point 
and the constant c has to be determined by the ap- 
propriate boundary conditions, m satisfies the rela- 
tions of a propagator, m(u,x,x') — m(u,x' jx)" 1 and 
m(it, x, x")m( ) — m(u, x, x'). Now we determine 

the constant c for a system of size [— d s ,d]. We as- 
sume specular reflection at two boundaries at x = — d s , d 
and ideal interfaces between different materials inside 
the system. The appropriate boundary conditions are 
f(v x ,Vy,x = —d s ,d) = f(—v x ,Vy,x — ~d s ,d) and conti- 
nuity at the internal interfaces. The same conditions are 
valid for a\ and a\. This leads to 



c(V X ,Vy 



2e^ 

v x 



m(v x , d, x') + m(—v x , d, x') 



ds m(v x , d, -d s ) - m(-v x , d, -d s ) 
x[l- go (v x ,x')}A(x')dx' . (20) 



The current is determined by Eq. ([}]), using the Green's 
functions (|l^) , expressed by the solution ( |l8| ) . We obtain 
the following general result for the linear current response 
in functional dependence of the vector potential, 



3v(x) 



K(x,x')A(x')dx' 



(21) 



where the kernel K(x, x') is given by 
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K(x,x') 



2 2 



7T 



w>0 



[1 + g (u,x)][l -g (u,x')] 



Q(x — x')m(u, x, x') + Q(x' — x)m(—u, x, x') (22) 



m(—u, x, d)m(u, d, x') 



m(u, x, —ds)m{—u, —d s , x') 



1 — m(u, d, —d s )m(—u, —d s , d) 1 — m(u, d, —d s )m(—u, —d s , d) 

m(—u, x, d)m(u, d, —d s )m(—u, —d s , x') m(u, x, —d s )m(—u 1 —d s , d)m(u, d, x') 



1 — m(u, d, —d s )m(—u, —d s , d) 



1 — m{u, d, —d s )m{—u, —d s , d) 



Equation (p2) gives the exact linear-response kernel of 
any quasi-one-dimensional system, consisting of a combi- 
nation of normal and superconducting layers extending 
from x = —d s to x = d. The kernel is expressed in 
terms of the quasi-classical Green's functions in absence 
of the fields, which may be specified for the particular 
problem of interest. We note two characteristic features 
of Eq. (H|): The factor (1 — g (u,x')) measures the de- 
viation of a quasi-classical trajectory from the normal 
state go = 1, which is inert to a magnetic field. The 
propagator m(u, x, x') shows up in six summands which 
represent all the ballistic paths from x to x', accounting 
for multiple reflection at the walls at — d s and d. Thus 
the first two summands connecting x and x' directly con- 
stitute the bulk contribution, while the additional four 
summands are specific to a finite system (assuming spec- 
ular reflection at the boundary). We note that a form 
similar to (]2^) may be derived for non-ideal interfaces 
between the normal and superconducting layers, if the 
appropriate boundary conditions following from Ref. |2^ 
are taken into account (these boundary conditions are 
only valid if the distance between two barriers is larger 
than the mean free path). 

For illustration we reproduce the current response 
of a half-infinite superconductor. Setting d = and 
d s — > oo, the solution of the Eilenberger equation (Q) 
takes the simple form go = us/VL, /q = /J = A/f2, where 
O = (A 2 +W 2 ) 1 / 2 . Inserting in (p2h we obtain the linear- 



response kernel 



Ks(x,x') 



2 pl v-^ A 2 f 1-u 2 / 



UJ>0 



(23) 



- ( 2o+i)M_£a +e(2 o+i)2±£: 



which describes the current response of ag-arbitrary su- 
perconductor, as first derived by GorkovEa, which here 
additionally includes the effect of the boundary. For 
fields varying rapidly spatially we arrive at,a non-local 
current-field relation of the Pippard-typeo, while for 
slowly varying fields the kernel can be_integrated out in 
Eq. producing the London resultEJ. We recall here 
certain generic features of this kernel, which are of impor- 
tance below. In a dirty superconductor (fi <C 1/r) the 
range is given by the mean free path I — vft. In a clean 
superconductor (1/r <c 0), the range is roughly given by 
the coherence length £o and is thus nearly temperature- 
independent. 



IV. MAGNETIC RESPONSE 

For the NS system we consider in this paper, see Fig. ||, 
the kernel ( [22] ) may be simplified using m(u, x, ~ oo) — * 
and m(—u, — oo, x) — > as d s — > oo (u > 0). The linear- 
response kernel takes the form, 



K(x,x') = -^TY] f du^^(l+g (u,x))(l-go(u,x')) 



Q(x — x )m(u, x, x ) + Q(x' — x)m(—u, x, x) + m(—u, x, d)m(u, d, x) 



(24) 



The magnetic response of the proximity system fol- 
lows from the self-consistent solution of Eq. (|2l]) and the 
Maxwell equation 

-^A(x) = -Anj v (x) . (25) 

As boundary condition we use -^A(x)\ x= d = H, where 
H is the applied magnetic field, and A(0) — 0, neglecting 
the penetration of the field into the superconductor. The 
inclusion of the field penetration into the superconductor 



leads to corrections ~ Xs/d to p, which is a small ratio 
for typical proximity systems. Here As is the effective 
penetration depth of the superconductor, including the 
nonlocal or impurity effects. The magnetic response of 
the normal-metal layer is measured by the screening frac- 
tion p = —Airx = 1 — A(d)/Hd, which gives the fraction 
of the normal-metal layer that is effectively field free. It 
is given by the susceptibility x, which is equal to the ra- 
tio of the average magnetization to the applied magnetic 
field. 



G 



The general properties of the kernel (|22|) are char- 
acterized by both the decay (range) of the propagator K(x, x ) = 
m(vF,x,x') and the amplitude of the prefactor (1 + 
go)(l — go) which determines the degree of non-locality of 
the relations (|2l]). The inverse decay length of the prop- 
agator is proportional to the off-diagonal part of the self- 
energy A and the prefactor is related to the superfluid 
density. We discuss below how, in the proximity effect, 
the range of the kernel varies from infinity to I and £t, 
exhibiting a strong temperature dependence, which leads 
to non-trivial screening properties. Furthermore, the su- 
perfluid density introduces an additional length scale in 
the problem: the London length Xn , which becomes cru- 3 ( x ) 
cial for the distinction of various regimes. 



2 2 

e • Pf 



(30) 



ui>0 
I*-*'! 





2d-z-z' 



+ e 



The kernel is factorized in a part containing the temper- 
ature dependence and a part which is responsible for the 
non-locality. The current is then expressed as 



K d {x,x')A{x')dx' . (31) 



Att\ 2 ( Xi T) 



A. Clean limit 

A special case is the clean normal metal (I —> oo). Here 
the range is infinite and the current-field relation is com- 
pletely non-local. It follows from (^2)) that it is necessary 
to have impurities in a normal metal to get a finite range 
of the kernel. In the limit d > (o the current may be 
written as 



1 



Jclcan 



4irX 2 (T)d 



A(x)dx 



(26) 



This defines a temperature-dependent penetration depth 
that can be given explicitly in the limits T — and 
T > T A : 



A 2 (T) 



Aire n e 



= : X 



N ) 



T = 



-e t a 



(27) 



12TA 

Solving Maxwell's equation we find 

A(x)=Hx(l- 3 -- d ( M -*> 



T>T A . 



A{?>X 2 {T)+d 2 )' 



and the screening fraction 



4+12A 2 (T)/d 2 ' 



(28) 



(29) 



In the limit X(T) <C d the screening fraction is 3/4, thus 
the screening is not perfectly diamagnetic. The magnetic 
field inside the normal metal is B{x)/ H = l—2p(l—x/d), 
showing the effect of overscreening for p > 1/2, where the 
field reverses sign inside the normal metal. 



B. Dirty limit 

Using the fact that the zeroth-order Green's function 
is nearly isotropic and varies on a scale £d(T) 3> I, we 
find for the kernel (M) 



The local penetration depth X(x,T) is defined as 
1 iwr 



(32) 



and the temperature-independent part of the kernel is 
given by 



K d (x,x) = 



,\x — x'L ,\x — je'1 . 
EtC—r-t) - EsC—r-t) (33) 

,1d — x — x' , 2d — x — x' , 

-£ X ( ] )-E 3 (- 



4/ 



I 



In this formula E n (z) = t~ n exp(—zt)dt is the ex- 
ponential integral. For X(x,T) ^> I the vector poten- 
tial may be taken out of the integral in Eq. (|3l]) and 
the spatial integral yields the well-known local pirrent- 
vector potential relation used in Usadel theoryEH We 
note that for A <C Co there may exist a region (see Fig. |^) 
where the Green's functions are nearly isotropic, and in 
absence of the field are given by Usadel theory, but the 
current response is nonlocal. To put limits on the validity 
of the local relation, we consider the approximate form 
(/) ~ exp[— x(2u/D) 1 ^ 2 ] to determine the local penetra- 
tion depth. As a result we find 



X(x,T) 



Ajvf if £r>(T)»d 



(34) 



To achieve locality we need to have I < X(x,T) in the re- 
gion, where the screening takes place. For T <C D/d 2 this 
means I <§C X(d) leading to the condition I 2 <C Xjyd. For 
T 3> D/d 2 screening takes place at x w £d and we have 
I 2 <C Xn£,d(T). The local penetration depth at the outer 
boundary can be small compared to d. In that case the 
screening fraction p = 1 — X(d)/d can reach practically 
unity. 

At high temperatures T S> -E/rh, the inverse penetra- 
tion depth is exponentially suppressed on a scale of the 
temperature-dependent coherence length £d(T). This 
length defines the screening region and consequently 
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p(T)oc 



d 



cxT' 



-1/2 



(35) 



This result has already .been obtained on the basis of 
Ginzburg-Landau_theoryQ and numerically confirmed us- 
ing Usadel theorycj. We expect that non-local screening, 
which may be taken into account using Eq. (|3^) , will only 
lead to quantitative corrections to Eq. ([55]). We will not 
consider this here but concentrate on the more interest- 
ing case in which non-locality gives rise to a qualitatively 
different picture. 



C. Arbitrary impurity concentration: numerical 
results 



As has been shown in the last two subsections, there 
are two main differences in the observable properties of 
the induced screening in the clean or the dirty limit. 
Firstly, the saturation value of p in the dirty limit can 
reach practically unity, whereas in the clean limit it is 
limited to 3/4. The analytic behavior at high tempera- 
tures is quite different too. In the dirty limit p shows an 
algebraic behavior oc T -1 / 2 , whereas in the clean limit 
we find p oc exp(— 2T/Ta). From a theoretical point of 
view these two limits are characterized by a completely 
nonlocal constitutive relation in the clean limit and a lo- 
cal relation in the dirty limit. In this section we will 
investigate the magnetic response in the regime between 
these two extreme cases. 

To calculate the diamagnetic response, we have evalu- 
ated the integral kernel in Eq. ( p4| ) numerically using the 
results from the Sec. || and solving Maxwell's equation 
by a finite-difference technique. Therefore, the parame- 
ters entering the calculation are I / d and Xn /d, assuming 
I > £o; an d ^r/d giving the temperature dependence. 

Magnetic field distributions for various impurity con- 
centrations and temperatures are shown in Fig. ^. Thick 
lines show the magnetic field and thin lines the current 
distribution inside the normal metal. Different graphs 
correspond to different mean free paths and the curves 
inside each graph to different temperatures. In all curves 
we have chosen Ajv = 0.003d. All these curves clearly 
show deviations from the clean-limit behavior, where the 
field decays linearly and the current density is spatially 
constant. Larger impurity concentrations make it possi- 
ble to localize the current on a length scale smaller than 
system size. Obviously, this scale is not given by the 
mean free path, but can be considerably smaller. Later 
we will show what determines this length scale. Whether 
the localization of the current increases or decreases the 
screening fraction depends on temperature. Note that 
for all parameters chosen the field is overscreened, which 
is the signature of a nonlocal constitutive relation. 
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FIG. 4. Magnetic induction (thick lines) and current den- 
sities (thin lines) in a proximity layer for different mean free 
paths. The different curves in each graph correspond to tem- 
peratures of T/Ta = 0.04 (solid line), 5 (short-dashed line), 
and 8 (long-dashed line). 
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FIG. 5. Numerical results for screening fraction of the 
normal metal layer for Ajv = 0.003d. The clean limit is indi- 
cated by a thin line reaching p = 0.75 for T — » 0. Even a very 
large mean free path of 10 4 d leads to an enhanced screening 
at low temperatures. For smaller mean free paths (but still 
> d) the screening is enhanced at high and low temperatures. 
At the smallest mean free path Q.ld the system is in the dif- 
fusive regime leading to a completely different temperature 
dependence. 



The screening fraction as a function of temperature is 
shown in Fig. |5| The different curves are for the clean 
limit and for mean free paths l/d = 10 4 , 10, 1, 0.1. We see 
that a finite impurity concentration has strong influence 
on the screening fraction, even if I > d. It can either in- 
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crease or decrease the diamagnetic screening, depending 
on temperature. 

For the interpretation of these results we first consider 
the case I = O.ld. The lower-right graph in Fig. |J shows 
that the screening is nearly local, since overscreening is 
rather small. The local screening strength depends on the 
local superfluid density. At low temperatures T <C Ta 
the superfluid density at x = d is finite and the field is 
screened exponentially, leading to a screening fraction of 
nearly unity. A higher temperature suppresses the super- 
fluid density and the field penetrates to the point where 
the density is large enough to screen effectively. On the 
other hand, the locality of the kernel allows the system to 
screen even if the superfluid density is suppressed nearly 
everywhere. The screening is then enhanced in compari- 
son to the clean limit. In Fig. ||this appears at T s=s 6Ta- 

Let us now consider a mean free path of order or much 
larger than the sample size. Even for / = 10 4 d we see a 
deviation from the clean-limit expression at low temper- 
atures. For / = lOd and / = d screening is enhanced in 
comparison to the clean limit at low and high temper- 
atures. Only in an intermediate regime, i. e. around 
T = 5Ta in our case, p is similar to the clean limit 
screening fraction. A qualitative understanding may be 
gained from looking at the constitutive relation in the 
limit / 3> d. In the limit TcTa the zeroth-order Green's 
functions are given by the clean-limit expressions fla) . We 
approximate the kernel ( |24| ) by 



K(x,x') 



1 



8ir\ 2 (T)d 



(36) 



Since I ^> d, the exponentials may be expanded to first 
order. As a result, we obtain two contributions to the 
current 



-1 



Jclcan 



Jimp(^) 



4ir\ 2 (T)d J 
1 

8nX 2 (T)d Jo 



A(x)dx (37) 

(38) 

d \x -x'\ + 2d-x-x' . . , , 
1 ■ — A(x')dx' . 



When will deviations from the clean limit become im- 
portant? It is clear that the im pur ities cannot be ne- 
glected, if ( ^8|) is comparable to (|37|) . We estimate this 
by calculating the two contributions to the current using 
the clean-limit vector potential (p8|). Comparing the two 
contributions, we find that impurities can be neglected, 
if 



Aeff(T) = ZJ\ 2 (T)l » d 



(39) 



This equation defines a new length scale, the effective 
penetration depth A e fT, which determines the validity of 
the clean-limit magnetic response. For the clean limit to 
be valid at T — the condition A e ff(0) > d has to be 
fulfilled, since in this case the screening takes place on 
the geometrical scale d. In the case A o ff(0) -C d the field 



is screened on a scale A e ff and the screening fraction is 
strongly enhanced in comparison to the clean limit. Nev- 
ertheless, the clean-limit behavior reappears at higher 
temperatures, since A e ff(T) grows with temperature. 

For T 3> Ta the deviations from the clean limit are re- 
lated to deviations of the zeroth-order Green's function 
from the clean-limit expression due to impurity scatter- 
ing. The correction to g, given in Eq. po|), leads to a 
finite superfluid density in the vicinity of the supercon- 
ductor via the factor 1 — g(x') in the kernel. The range 
of the propagator is modified by the correction to /+, 
leading to 



m(x, x ) 



exp (i 



dx 



• {fix")) 
p{x") 



exp 



(40) 

Thus, the range of the kernel is now given by which is 
strongly temperature dependent. Summarizing, we find 
for the current 



j(x) 



-1 



Aeff(0) S 



dx e 



It 



A(x') , 



(41) 



again showing the importance of the new length scale A c ff . 
In the limit A c _t(0) 3> £r the field cannot be effectively 
screened on the scale £t , leading to a vanishing screening 
fraction. If A o ff(0) <§C £t the field can be screened on a 
length scale smaller than £t and the screening fraction 
will be finite. 

It is therefore evident that the interplay between lo- 
cal and nonlocal physics is of crucial importance for the 
screening behavior of a normal-metal proximity layer. 
The most interesting regime occurs for I > d, where a 
transition between different screening behaviors may be 
observed by varying the temperature, see Fig. ^[ 

We note that the screening fraction is a non-monotonic 
function of the mean free path. At low temperature, with 
increasing mean free path (i.e. increasing purity), the 
screening fraction is reduced rather than enhanced. As- 
suming a temperature-dependent scattering mechanism 
with decreasing mean free path as a function of temper- 
ature, such as electron-electron or electron-phonon inter- 
action, we might speculate to observe a non-monotonic 
(i.e. re-entrant) behavior of the susceptibility (here the 
smallness of the scattering rate is compensated by the 
high sensibility of the non-local current-field relation). 
However, as is evident from Eq. (|l9|), the largest off- 
diagonal self-energy (A) which includes e.g. impurity 
scattering will provide a (low-temperature) cutoff for this 
behavior. 

Finally we comment on the effect of a rough bound- 
ary. For T Ta the Green's functions are independent 
of the boundary condition at x — d. In this case a finite 
screening fraction can only be due to impurity scatter- 
ing, see Fig. |. For T <T A the screening behavior will 
be strongly affected by a rough boundary. This makes 
it possible to distinguish between clean samples with a 
rough boundary and samples containing impurities. 







V. CONCLUSIONS 

We have investigated the diamagnetic response of a 
proximity layer for arbitrary impurity concentration us- 
ing the quasiclassical theory of superconductivity. We 
found a variety of different regimes in which the physics 
is different from the previously studied clean and dirty 
limits. 

We have first investigated the proximity effect in the 
absence of fields, distinguishing three different regimes, 
see Fig. ||. In the ballistic regime, the validity of the 
clean-limit solution is restricted to I ^> dcxp(2d/^T) for 
£t *C d and to I ^S> d for £r 3> d. The last condition 
is the consequence of the suppression of the density of 
states for to <C v-p/d, which enhances the effective mean 
free path to ~ Ivp/ud. In the diffusive regime, we found 
that the validity of the Usadel equation (dirty limit) de- 
pends on the superconductor as well as the normal metal, 
and is thus restricted to I <C £o and Icd,^. The first 
condition is due to the fact that the induced supercon- 
ducting correlations are strongly anisotropic for a clean 
superconductor (I 3> £o)> even if the motion is diffusive 
in the normal metal. The intermediate diffusive regime 
(Co -C I <C d) is not covered by these two cases. Here the 
full Eilenberger equation has to be solved, which requires 
a numerical analysis, see Fig. |[ 

To study the magnetic response of the proximity layer 
we have derived explicit expressions for the general 
linear-response kernel ( p2| ) for an NS sandwich. This 
derivation may easily be generalized to systems such as 
Josephson junctions or unconventional superconductors. 
We have used this linear-response kernel to study the 
magnetic response of the proximity system at arbitrary 
impurity concentrations. The nonlocal current-field re- 
lation is shown to have non-trivial consequences on the 
screening behavior of the normal metal. In the ballis- 
tic case, we found the clean-limit theory to be restricted 
further by d < {\ 2 {T)l) 1 /' i = X c g, A e ff giving the pen- 
etration depth for the non-local current-field relation. 
If A c ff > d, the screening takes place on the geometric 
length scale d, leading to a saturation at the screening 
fraction of 3/4 at low temperatures. If A e s < d, the finite 
(even though large) mean free path strongly enhances 
the screening. Thus for typical samples with A at <C d 
even a mean free path I 3> d cannot be neglected, i.e., 
the clean-limit behavior is practically unobservable. At 
large temperatures T 3> Ta, a finite impurity concen- 
tration reduces the range of the linear-response kernel 
to £t, again enhancing the screening. Furthermore, the 
screening fraction may serve to distinguish between sam- 
ples with bulk impurities rather than a rough boundary, 
since a nonzero screening fraction at large temperatures 
is only due to bulk impurity scattering. In dirty systems, 
where the zeroth-order Green's function is well described 
by the Usadel approximation, the current-field relation 
can still be non-local. We have shown that the appli- 
cability of the local current-field relation is restricted to 
I 2 -C \ N d for T < E Th and I 3 < X 2 N f D {T) for T > E Th . 



This shows that in the presence of magnetic fields some 
caution is needed in applying the Usadel theory. 
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